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The Ford-Kac-Mazur formalism is used to study quantum transport in (1) electronic and (2) 
harmonic oscillator systems connected to general reservoirs. It is shown that for non-interacting 
systems the method is easy to implement and is used to obtain many exact results on electrical and 
thermal transport in one-dimensional disordered wires. Some of these have earlier been obtained 
using nonequilibrium Green function methods. We examine the role that reservoirs and contacts 
can have on determining the transport properties of a wire and find several interesting effects. 
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I. INTRODUCTION 

There is considerable current interest in the prob- 
lem of transport through various nanoscale devices both 
from the fundamental and from applied points of view. 
In this connection, Kubo's transport formulas have to 
a large extent been superseded by different formalisms 
in the spirit of Bardeen's tunneling model^. The Lan- 
dauer formula^ (LP) and the Keldysh technique'^, quan- 
tum Langevin equations*, C* algebraic formulas^ and 
generalized scattering theory ideas® have been developed, 
allowing one to study systems in steady state arbitrar- 
ily far from the linear region where Kubo is applicable. 
There is also considerable experimental activity involv- 
ing resistive elements, such as quantum dots, STM tips, 
single walled nanotubes and insulating nanowires, often 
coming up with unexpected physics^"^. 

The most popular alternative to Kubos formulas is the 
LF, proposed in 1957^. Since then several derivations 
of the LF have been given^" and this has led to a good 
understanding of the formula. A large number of experi- 
ments are interpreted successfully on the basis of LF. The 
quantum of conductance /h has been understood as a 
contact resistance which arises due to the squeezing of 
the reservoir degrees of freedom into a single channcl^^'^^. 
While a physically careful statement of the conditions for 
validity of the LF can be found in Ref(^^), we believe that 
a detailed mathematical theory of the role of reservoirs 
and the nature of the coupling between the wires and 
reservoirs does not exist. The role of the idealized reser- 
voirs has been to serve as perfect sources and sinks of 
thermal electrons. This clearly will not be satisfied in all 
experimental conditions and it is necessary to have a bet- 
ter microscopic understanding of reservoirs and contacts. 
There has been some work'^'^'®'^^d4 ^j^jg direction but, 
to our knowledge, a detailed understanding of the role of 
reservoirs is still lacking. 

In this paper we adapt a formalism that was developed 
by Ford, Kac, and Mazur-'^^ (FKM) and model reser- 
voirs as infinite non-interacting systems. This method 
was originally devised to study Brownian motion in cou- 
pled oscillators"'^^ and was later extended to a general 



study of the problem of a quantum particle coupled to a 
quantum mechanical heat bath^® . In this approach reser- 
voirs are modelled by a collection of oscillators which are 
initially in equilibrium. The reservoir degrees of free- 
dom are then eliminated leading to quantum Langevin 
equations for the remaining degrees of freedom (the sys- 
tem). Thus the reservoirs can be viewed as providing 
sources of noise and dissipation into the system. The 
FKM formalism is thus very direct to interpret and, as 
we shall demonstrate, is more straightforward to apply 
than other methods of treating open quantum systems 
such as the Caldeira-Leggett^^, Keldysh^* and scattering 
theory®. Quantum Langevin equations have earlier been 
used in the context of transport in mesoscopic systems 
and have helped in the understanding of some experimen- 
tal data*'^^. The FKM approach was also used earlier 
by O'Connor and Lebowitz^" in studying classical heat 
transport in disordered harmonic chains and our anal- 
ysis here closely follows theirs. Here we use the FKM 
approach to make a detailed study of quantum trans- 
port in disordered electronic and phononic systems. For 
very general reservoirs we obtain exact formal expres- 
sions for currents and local densities in the nonequilib- 
rium steady state. We find that for a special type of 
reservoir, the ideal Landauer result (where the conduc- 
tance is expressed in terms of the transmission coefficient 
of one-dimensional plane waves) follows exactly, while for 
general reservoirs they need to be modified. We examine 
in some detail the effect on transport properties that the 
choice of reservoirs can have and find a number of inter- 
esting effects. For example in the electron case we find 
that imperfect contacts can lead to an enhancement of 
conductivity. In the phonon case we find the surprising 
result, earlier noted for classical systems, that the heat 
current J in a long disordered wire decays with system 
size N as J ^ where a depends on the low fre- 

quency spectral properties of the reservoirs. 

The paper is organized as follows: In Sec. I we present 
the formalism and results for transport in the one di- 
mensional Anderson model. In Sec. II we present the for- 
malism and results for transport in disordered harmonic 
chains. We end with a discussion in Sec. III. 
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FIG. 1: A disordered system connected through ID leads to 
reservoirs at different potentials and temperatures. 



II. TRANSPORT IN THE ONE-DIMENSIONAL 
ANDERSON MODEL 

A. The formalism and main results 

The set-up: We wish to study conduction in a disor- 
dered fermionic system connected to heat and particle 
reservoirs through ideal ID leads [see Fig. 1]. We con- 
sider a tight-binding model and for simplicity we take 
the system and leads to bo ID while the reservoirs are 
quite general. We use the following notation: the indices 
l,m denote points on the system or leads, greek indices 
XjP or X',fi' denote points on the left or right reservoirs 
respectively, finally p, q denote points anywhere. Thus 
c; (/ = 1, 2...N) denotes lattice fermionic operators on the 
(system -|- lead), c\, cy (A, A' = 1,2...M) denotes oper- 
ators on the left and right reservoirs. The CpS' satisfy the 
usual anticommutation relations {cp,Cq} = 0; {cj,,cj} = 
0; {4, Cq} = Spq. Out of the N = Ns + 2Ni, sites the first 
and last Ni sites refer to the leads while the middle A^^ 
sites refer to the system. The Hamiltonian for the entire 
system is given hy H = + V + Vint, where 



of the system and leads is given by (for t > t): 

C\ = ici — iv\C\ -|- i^Ca 
ci = ?:(q_i +Q+i) - W;Q (2<Z<iV-l) 
CAT = iCE-\-ivj^Cf^ ^iica'. (1) 

The equations at the boundary sites involve reservoir op- 
erators, Ca, Ca'. Using the equations of motion of the 

reservoir variables wc can replace these reservoir opera- 
tors by Langevin type terms. The equations for the left 
reservoir are given by (for t > t): 



c\ = -iTxuC„ (A 7^ a) 



(2) 



N-l 



N 



H° = -Y^ (4ci+l + 4_^^Cl) + Vic\ci +Y^fxJyC^ 



l=\ 



l=\ 



The first part of refers to the system and leads, while 
T and T' describe the two reservoirs. The contact be- 
tween the reservoirs and leads is given by the intercon- 
nection part V . The interaction part Vint can be added 
perturbatively, and we return to its inclusion later in the 
paper. We will consider a system with onsite disorder 
and so choose the onsite energies w;, I — N1 + I...N1 + Ng, 
from some random distribution. At sites belonging to 
the leads [/ = 1, 2, ...Ni, Ni + Ns + 1, ...A^], assumed to be 
perfect conductors, we set vi = 0. At some time t < t 
in the remote past, the two reservoirs are isolated and in 
equilibrium at chemical potentials /z and /z' and inverse 
temperatures (3 and /?' respectively. At t = t, we con- 
nect the reservoirs to the two leads and evolve the system 
with the Hamiltonian H. We study the properties of the 
nonequilibrium steady state, reached after a long time. 
The Heisenberg equations of motion for the operators 



This is a linear set of equations with an inhomogeneous 
part given by the term i'yci and has the general solution 

POO 

exit) =iY^9Ut - r)cAr) - / dt'g+Jt - t'){jc,{t')) 
where gl^t) = -i0{t)Y,^n{X)rn{^)e-''-' . 

n 

Here Vn(A) is the single particle eigenstate of the left 
reservoir, with energy e„, and n runs over all states. 
We need Ca{t) which we note has two parts. The first, 
^i-it) = ^X^iyfl'a!/(* ~ ''")ci/(t), is like a noise term whose 
statistics is determined by the initial conditions of the 
reservoir. Initially the reservoirs are in thermal equi- 
librium and the normal modes c„ = X);^CAV'n(A) sat- 
isfy (cJj(r)c„/(T)) = (5„„//(g„,/i,/3), where / is the Fermi 
distribution / = l/[ef>i'^^^i'') + 1] and (6) = Tr[dp] 
where p is the reservoir density matrix at time r and 
Tr denotes a trace over reservoir variables. The sec- 
ond part of Ca{t), —7 dt'g^^{t — t')ci{t'), is dis- 
sipative in nature. Defining the Fourier transforms 

and h{uj) = ^ dth{t)e''^^, and taking limits M —>■ 00 
and T — > —00 we get: 

Ca{uj) = h{u) - 7ff+„(a;)ci(a;) 
{h\u))h{uj')) = I{uj)5{uj - J), 
I{uj) = pa{uj)f{uj), 



- Pa 



(3) 



where Pa = Ylm I '^n{ct) P (5(w — e„) is the density of 
states at site a. The third equation above is a statement 
of the fluctuation-dissipation theorem. Similarly for the 
right reservoir we get c^'ioj) = h'{u) - j'g'^,^,{uj)cN{uj), 
with the noise statistics of h'{cu) determined by /x' and /3' . 
Also h and h' are independent so that {h"^ {Lo)h' {lo')) = 0. 
We now Fourier transform the system equations and plug 
in the forms of €0(0;) and Ca'(a;) to get the following 
particular solution: 



ci 



J — < 



(LjZi^{uj)hm{uj)e 



-iujt 



(4) 
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Zlm = ^Jm + Mm 

= -<J/,m+l - + ("i - 

/l/ = 7ft.(w)(5/,i +7'/l'(w)(5;,Ar. 

With this formal solution and the known properties of 
the spectral functions /i(a;), h'{uj), gaai^) g^'a'i^)' 
we can now compute various physical quantities of inter- 
est. Specifically we shall be interested in the electrical 
and thermal currents and the local particle and energy 
densities. The operators corresponding to particle and 
energy densities are given by 

hi = c\ci 

ui = -{cjci+i + cJ^iQ) + ^{vicjci + vi+icj^-^ci+iib) 

while the corresponding current operators j" and 
are defined through the conservation equations dn/dt + 
dr/dx = and du/dt + dj^'/dx = 0. We get: 

V t t \ 

jT = -^{cl,Cl - c|q+2) + '-^Or+i +.57) (6) 

We now calculate the steady state averages of these 
four quantities. We introduce some notation and state 
a few mathematical identities. We denote by Yi^m the 
determinant of the submatrix of Z beginning with the 
Zth row and column and ending with the mth row and 
column. Similarly D;^m denotes determinant of the sub- 
matrix formed from $. The following results can be 
proved: («) Yi,n £'i,w+7^ff^a^2,JV+7'^5^'a'-'^i,JV-i + 
7V'5+„,9i.'^2,Ar-i; (ii) Z-^' = Zf,' = 

5^i+l,Ar/il,JV; (iii) Dl,n-lD2,n - ■Dl,„£'2,n-1 = 1- 

Particle and Heat Currents: The expectation value of 
the current operators, using Eqns. (3,21), gives: 

/oo ^ 
du;Im[ Zf_^\:MZ^ri^)Ir{o^)] 

-°° r=l,JV 

/oo 
d^imi z^+\:A^)z^ri^)ir{u^)m 

where Ii = 7^/; In = Y'^I'- In the case of the heat 
current we take I to be on the leads so that vi = 0. 
Using the various identities stated earlier we can show, 
as expected, that these are independent of I and reduce 
to the simpler expressions 

/oo 
du;J{co)[f{u;)-nu;)] 
-OC 

/OO 
dcoujJ{co)[f{co)-f'iuj)] where 
-00 

J{uj) = 27r7V'pa(a;)Pa'(a;)/|Fi,;v p. 

These can be expressed in terms of the retarded Green's 
function G+(a;) = {ui + ie — H)~^. This satisfies G+ = 



g++g+VG+ where 5+ = (a;-|-ie-iJ°)~^ These can be 

solved to give : 

G^m = [9trm-^'^9tJ9ti9Nm- 9Nigtm)]/Z whcrc 

'7 1 2++ /2+ + 

Z = l-"f gugZa - 7 9nn9Z'c' 

Wi^9aa9l'a'{9tl9tlN " ff]JV.9^l) 

Let gim = i?e[(7;'^] denote the real part of the system's 
Green function. It is easy to see that giN = gni = 

— ^/Di^N', gn = —D2,N / Di^N'-i gNN = —Dl,N-l/Dl,N- 

Using these and the Jacobi identity gngNN — giNgNi = 
D2,N-i/Di^N we get 1/| Yi^n P = Gn^j^Gn]^^ where 
Gn+ is a modified Green function obtained from G+ by 
replacing all system green functions by their real part. 
We then get the particle current in a form similar to those 
obtained by Meir and Wingreen^^ using the Keldysh for- 
malism and by Todorov et al^ using time-independent 
scattering theory. Their results differ from ours in that 
they are expressed in terms of G+ instead of Gn+. The 
case of insulating wires treated by Caroli et aP also fol- 
lows from our results. 

Scattering states: It is instructive to write the currents 
and densities in terms of properties of the single-particle 
scattering states of the full Hamiltonian H ( possible 
when interactions are absent). Let iIj^^{ui) and ip^^{u)) 
denote the jth unperturbed wave functions with energy 
CO, of the left and right reservoirs respectively. Let a^^ 
and o^^ denote the amplitude at site p of the jth right 
and left moving states obtained by evolving the unper- 
turbed levels with the full Hamiltonian. We then get: 
af^ = Kf.^'y^iHw); = Kf^Wi^l'^i^). The cur- 
rents and densities arc given by j" = i{af^-^^ai — a*a;+i), 
Hi = a^ai etc. Using these we find that J{ijj) is simply 
the total transmitted current for all waves with energy 
w. Also the particle density is given by 

/oo 
da;[pf'(a;)/(a;)+pf(a;)/'(a;)], 
-00 

where Pi = la/^P is the total particle density at a 
point I due to all right moving waves with energy lo and 
is due to left movers. Note that the currents and 
densities do not have the simple Landauer form since 
J{uj) depends not only on the system but also on bath 
and contact properties. The spectral properties of the 
baths enter into the expressions in a nontrivial way and 
one cannot separate the contributions of the system and 
the baths. 



B. Ideal reservoirs and contacts:The Landauer 
case 

This corresponds to the case where 7 = 7' = 1 
and the reservoirs themselves are semi-infinite exten- 
sions of the one dimensional leads. This results in re- 
flectionless contacts between the reservoirs and leads. 
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The reservoir wave functions and energy eigenvalues are: 
ipn{X) = [2/{M + l)Y/^sin{kX), e„ = -2cos(fc) where 
k = mT/{M + 1) with n = 1,2,...M. The leads are 
connected at the end of the reservoir chains so that 
a = a' = 1. We then get the following reservoir spectral 
functions: 

TT 

I{u}) = Iwl > 2 

5i"aM = -e''', -2<a; = -2cos(A;) <2 

= u;/2- sgn{uj){Lj'^ /A- 1)^/'^, > 2 (8) 

Wc have similar expressions for the right reservoir. Let 
us use the notation that if sites iV; + I and Ni + m be- 
long to the system then we write Yni+i^Ni+tti = yi,m, 
DNi+i,Ni+m = di,™.- It can be shown that the transmis- 
sion probability of a wave with momentum k across the 
system is given by 



T=i£^ where 



(9) 



\yi,NA = \di,N, - e"'{d2,N, + di^N.,-i) + e'^''d2,NlM 

Note that in this case the transmission factor does not 
involve properties of the reservoirs and contacts. Also 
transmission is only by propagating modes which can 
be labelled by a real wave vector k (In general, non- 
propagating modes would also carry current and we 
would have integrate over all frequencies). Wc then get 
the following forms for the particle and energy currents: 

OT) = ^ [ dku{k)T{k)[f - f] 

Or) = ^ r dkv{k)e{k)T{k)[} - /'], where(ll) 

271- 7o 

v{k) = de{k)/dk = 2sin(fc). 

which are precisely of the Landauer form. 

In order to get the four-probe result we need to find 
the actual potential and temperature differences across 
the system. We imagine doing this by putting poten- 
tiometers and thermometers at points on the leads [A 
and B in Fig. (1)]. These measure the local particle and 
energy density on the leads from which one can compute 
the chemical potential and temperature. Wo note that we 
do not expect local thermal equilibration in this nonin- 
teracting system and so these are only effective potentials 
and temperatures. 

We start with the general expressions for densities 
[similar to Eqns. (7)] and after using the various determi- 
nantal identities we get ( for points / located on the left 
lead) an integrand which contains a factor sin^[fc(iVi — I)]. 
Assuming that Ni is large and / is not too close to the 
point of contact with reservoirs this factor can be re- 
placed by 1/2. We then get for the particle and energy 
densities: 



(ui) = ^ r dke{k){[2 - Tik)]f + T{k)f'}. (12) 

27r Jo 

We get similar expressions for densities at points on the 
right lead. The expressions in Eqns. (11,12) are identical 
to those obtained from semiclassical arguments, are true 
for ideal contacts, and lead to the usual four-probe for- 
mulas. The results of Eqns. (11,12) have been obtained 
earlier by Tasaki^ using the theory of C*— algebra. They 
can be easily extended to the case where the leads are 
still one-dimensional but the system is of more general 
form. Thus let the system consist of Ns points of which 
1 and Ng are connected to the two leads. Let us specify 
the system by the matrix (f> such that (fm ^ vi — lo and 
4'im = — 1 whenever two distinct points I and m are con- 
nected by a hopping element. Then all the above formu- 
las Eqns. (11,12) for currents and densities hold provided 
we evaluate them within the leads and use the appropri- 
ate expression for the transmission coefficient, namely 



T = 



4sin(fc)2F2 



\di,N, - e''=(d2,jv, + c?i,jv,-i) + e'2fed2,jV3-iP 



(13) 



{ni) 



■T{k)]f+T{k)r) 



where F denotes the determinant of the subniatrix 
formed from (f) by deleting the 1st row and iV^th column 
while d is, as before, but now constructed from (j). 



C. An application 

As an application we show how the experimental re- 
sults of Kong et al^ can be understood qualitatively using 
our results by assuming imperfect contacts. 

We consider again semi-infinite ideal reservoirs but 
make the contacts non-ideal by setting 7 = 7' = 0.9. 
As system we take a wire with a single impurity at site s 
(Thus Vs 7^ 0). The linear response conductance is then 
given by G = da;J(a;)/(w)[l — /(w)]- We evaluate 
this numerically at different temperatures for A'' = 100, 
s = 10 and Vs = 0.2 [Fig. (2)]. We see the following fea- 
tures: (a) a rapid oscillation of the conductance due to 
resonances with standing waves in the wire, (b) a slower 
oscillation due to standing waves formed between bound- 
ary and impurity and (c) washing away of the oscillations 
with increasing temperature. These features are qualita- 
tively the same as seen in the experiments in^. The over- 
all decrease in conductance with increasing temperature 
is presumably due to scattering by phonons and hence is 
not seen here. We have also plotted in Fig. (2) the con- 
ductance as given by the usual LF. Note that this does 
not give the oscillatory features. Thus imperfect contacts 
cannot be treated as resistances in series with the system. 
Another rather remarkable effect we see is the enhance- 
ment of the conductance as a result of introduction of 
imperfect contacts. Infact we can sec in Fig. (2) that at 
certain values of /U the conductance almost attains the 
ideal value l/(27r). Similar features are also obtained if 
we make the contacts ideal but take other forms of reser- 
voirs (e.g rings or two dimensional baths). 
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FIG. 2: Plot of conductance-versus-Fermi level, at three dif- 
ferent temperatures (T) , of a wire with a single impurity and 

imperfect contacts. The corresponding plot for perfect con- 
tacts at T = is also shown. The horizontal line is the ideal 
conductance Go = l/(27r) [in units of /%]. 



D. Interacting systems 

For this case the present approach readily yields to 
a pcrturbative treatment. For illustration consider the 
case where the Hamiltonian of the system (and lead) Hsl 
contains an interacting part and is given by 

Af-l N N-l 

Hsl = - ^ + cJ+iQ) + X^^'(c|q + A ^ riini+i, 

1=1 1=1 1=1 

while the reservoirs are still taken to be non-interacting. 
In this case, Eqns. 1 take the form 

ci = ic2 — ivici — iAn2Ci + i'yca 

ci = i{ci-i + ci+i) - ivici - zA(n/_i -|- n;+i)c;, 

2 < / < TV- 1 
Cjv = icN-1 — ivNCN — iAnjv-iCjv +ij'ca' (14) 

and, being nonlinear, can no longer be solved exactly. 
However it is straightforward to obtain a pcrturbative so- 
lution which, schematically, has the form c{uj) = Z^^h — 
AZ-i / dw' J duj" Z-^hZ-^hZ-^h + 0{lS?). The opera- 
tors for particle density and particle current remain un- 
changed and we can obtain their expectation values as 
a perturbation series using this solution. Another possi- 
bility would be to solve Eq. (14) using a self-consistent 
mean field theory. 



III. HEAT TRANSPORT IN OSCILLATOR 
CHAINS 

We now use the FKM method to study heat conduc- 
tion in quantum disordered harmonic chains connected to 
general heat reservoirs which are modeled as infinite col- 
lection of oscillators. There has been some earlier work 
on quantum wires^^'^^ which follow a similar approach 
but we give a more clear and complete picture and make 
some interesting predictions for experiments. 

As in the electronic case we obtain formal exact ex- 
pressions for the thermal current and show that, for a 
special case, they reduce to Landauer-like forms. We 
also analyse the asymptotic system size dependence of 
the current and show that, depending on the reservoirs, 
a long wire can behave either like an insulator or a su- 
perconductor. Our results should be useful in interpret- 
ing recent experiments^ on heat transport in insulating 
nanowires and nanotubes. They are also of interest in the 
context of the question of validity of Fourier's law in one- 
dimensional systems, a problem that has received much 
attention recently^^. A large amount of work on classical 
Hamiltonian systems seem to indicate that Fourier's law 
is not valid in one-dimensionsal momentum conserving 
systems. Our work here shows that this is true even in 
quantum-mechanical systems. 

A. Formalism and main results 

We consider a mass disordered harmonic chain con- 
taining N particles with the following Hamiltonian: 

where {xi} and {pi} arc the displacement and momen- 
tum operators of the particles and {m;} are the random 
masses. Sites 1 and N are connected to two heat reser- 
voirs (L and R) which we now specify. We model each 
reservoir by a collection of M oscillators. Thus the left 
reservoir has the following Hamiltonian: 

Hl = ^^ + Y.-Kirr,XiXm (16) 

1=1 l.rn 
M 52 2 M 

s=l s=l 

where Kim is a general symmetric matrix for the spring 
couplings, {Xi,Pi} are the bath operators and {Xi,Pi} 
are the corresponding normal mode operators. They are 
related by the transformation Xi = UyXg where 
17; s, chosen to be real, satisfies the eigenvalue equa- 
tion J2iKniUis = iOgUns for s = 1, 2...M. The an- 
nihilation and creation operators ag, a| are given by 
flg = (Pg — icjsXs)/(2a;s)^/^, etc. and = a|as is the 
number operator. 
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The two reservoirs are initially in thermal equilibrium 
at temperatures Tj^ and Tfj. At time t = t the system, 
which is in an arbitrary initial state is connected to the 
reservoirs. We consider the case where site 1 on the sys- 
tem is connected to Xp on the left reservoir while A'' is 
connected to Xp> on the right reservoir. Thereafter the 
whole system evolves through the combined Hamiltonian: 

Ht = H + Hl-\- Hr - kxiXp - k'xNXp,. 

The Heisenberg equations of motion of the system vari- 
ables are the following (for t > t): 

miXi = —[2xi~X2\-\-kXp 
mixi = -{-xi-i + 2x1 - xi+i) 1<1<N 
TTiNXN = -[-XN-i + 2xn] + k'Xp' . (17) 

We note that they involve the bath variables Xpy . How- 
ever these can be eliminated and replaced by effective 
noise and dissipative terms, by using the equations of 
motion of the bath variables. Consider the equation of 
motion of the left bath variables. They have the form: 



Xn = -KniXi n^p 
Xp = -KpiXi + kxi 



(18) 



This is a linear inhomogeneous set of equations with the 
solution 

Xn = J2[F,,i{t - r)Xi{T) + Gr,l{t - t)Xi{t) 

I 

dt'G„p{t-t')kxi{t') where (19) 

Fnl{t) = 6{t)Y,UnsUls 



Gnlit) = e{t)Y,UnsUl 



sin(<x's<) 



Thus we find that Xp (say) appearing in Eq. (17) has 
the form Xp{t) = h{t) + k dt'Gpp{t - t')xi{t'). The 

first part, given by h{t) = Ylii^piit - 't)^^^) + Gpi{t - 
t)Xi[t), is like a noise term while the second part is like 
dissipation. The noise statistics is easily obtained using 
the fact that at time t = t the bath is in thermal equi- 
librium and the normal modes satisfy {al{T)as'{T)) = 
f{uJs,PL)Sss'- Here / = l/(e^" — 1) is the equilibrium 
phonon distribution and (O) = Tr[pO] where p is the 
reservoir density matrix and Tr is over the reservoir 
degrees of freedom. We define the Fourier transforms: 
xiiu;) = J-J^^dtxiit)e^'^\ G+picj) = j:^^dtGpp{t)e^-' 

and h{uj) = ^ J^^dth{t)e"^*. Taking limits M ^ oo 
and T — > —00 we get: 



feG+ (w)a;i(w) 



Xp{uj) = hiuj, , ..^pp 
{h{ijj)h{uj'))=I{ijj)6{uj+J 
7(.) = »M 



- ib{(j) where (20) 



^i^) = E"2^f'^^'^ ~ ^(w-l-Ws)] 



Similarly for the right reservoir we get Xpi = /i'(w) + 
fc'Gp,p, (a;)xjv(w), the noise statistics of h'{u}) being now 
determined by /?'. The left and right reservoirs are inde- 
pendent so that {h{io)h' {io')) = 0. We can now obtain the 
particular solution of Eq. (17) by taking Fourier trans- 
forms and plugging in the forms of h{Lo) and h'{u)). We 
then get: 



/OO 
-OO 



iujt 



(21) 



Z = 4>lm- ^Im with 
$lm = -{^l,m+l + Sl,m-l) + - mi^'^)Sl,m 
Aim = Sl,m[k^G^{uj)Sl,i+k'^G+p,{w)5l,N] 

hi{(j) = kh{uj)di^i + k'h'{u>)di^N- 

We can now proceed to calculate steady state values of 
observables of interest such as the heat current and the 
temperature profile. We first need to find the appropri- 
ate operators corresponding to these. To find the cur- 
rent operator j we first define the local energy density 

2 2 ' 

^ + + ^{xi - xi+i)'^ . Using the current con- 



ui 



servation equation du/dt + dj/dx = and the equations 
of motion we then find that ji = [xixi^i + a;;_ii;)/2. 
The steady state current can now be computed by using 
the explicit solution in Eq. (21). We get 

Cji)= r d^{iu;)[k'Z^l{u)Z^_\^,{-u)I{u) 

J — OO 

+k'^Zl-^'^{^)Z^\^{-^)I'{uj)] (22) 

The matrix Z is tridiagonal and using some of its special 
properties (see sec.IA) we can reduce the current expres- 
sion to the following simple form: 



-L 



du) 



ojb{oj)b'{u}) 



d^J{u:){f-n 



if-f) 



(23) 



where J{uj) ~ k'^k''^ujb{Lu)b'{iu)/'K \ y^jv P has the phys- 
ical interpretation as the total heat current in the wire 
due to all right-moving (or left-moving) scattering states 
of the full Hamiltonian (system -I- reservoirs) . Such scat- 
tering states can be obtained by evolving initial unper- 
turbed states of the reservoirs with the full Hamiltonian 
(see end of scc.IA). As before we have denoted by Yi^m 
the determinant of the submatrix of Z beginning with 
the ^th row and c;olumn and ending with the mth row 
and column. Similarly let Di ^ denote the determinant 
of the submatrix formed from 
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B. Ideal reservoirs and contacts:The Landauer case 

For the special case when the reservoirs are also one di- 
mensional chains with nearest neighbor spring constants 
= 1 and the coupling constants k, k' are set to unity, 
we have: G+, = 0^,^, = e"*'' where w = 2sin(fc/2), 
I{uj) = /(w)sin(fc)/7r for \ uj \< 2 and /(w) = for 
I w |> 2. In this case Eq. (23) simplifies further and has 
an interpretation in terms of transmission coefficients of 
plane waves across the disordered system. We get: 

J = — J ^dwu;\tN{u;)\''if-f') where (24) 

4sin^(fc) 



\tN{uj)\ 



is the transmission coefficient at frequency uj. We have 
thus obtained the Landauer formula^ for phononic trans- 
port. It is only in this special case of a one-dimensional 
reservoir and perfect contacts that we get the Landauer 
formula. The reason is that, only in this case is the 
transmission through the contacts perfect, and this re- 
quirement is one of the crucial assumptions in the Lan- 
dauer derivation. Note that in Eq. (24) (i)the transmis- 
sion coefficient does not depend on bath properties and 
(ii) transmission is only through propagating modes. For 
general reservoirs where we need to use Eq. (23) the fac- 
tor J{u)) involves not just the properties of the wire but 
also the details of the spectral functions of the reservoirs. 
Thus the conductivity of a sample can show rather re- 
markable dependence on reservoir properties as we shall 
see below. The above Landaucr-likc-formula has earlier 
been stated in^^ and derived more systematically in^''. 
We note that in the high temperature limit T ,T' ^ oo 
Eq. (24) reduces to the classical limit obtained exactly 

:„20-22 



C. Asymptotic system-size dependences 

In the case of electrical conduction the conductance 
of a long disordered chain decays exponentially with sys- 
tem size as a result of localization of states. In the case of 
phonons the long wavelength modes are not localized and 
can carry current. This leads to power-law dependences 
of the current on system size as has been found earlier 
in the context of heat conduction in classical oscillator 
chains. A surprising result is that the conductivity of 
such disordered chains depend not just on the properties 
of the chain itself but also on those of the reservoirs to 
which it is connected. It can be shown^^ that the asymp- 
totic properties of the integral in Eq. (23) depend on the 
low frequency (w ^ l/TV^/^) properties of the integrand. 
This means that we will get the same behaviour as in the 
classical case. We summarize some of the main results: 

(i) The classical case where the reservoirs are them- 
selves one-dimensional. In this case we put k = k' = 1 



UJ = 2sin(fc/2). This was treated by Rubin and Grcer^^ 
and it was found that that J ~ j^72- Thus the ideal 
Landauer case will also show this behaviour. 

(ii) The case of reservoirs which give delta-correlated 
Langevin noise corresponds to taking k = k' = 1 and 
G+, = = —i-ycj. The classical case was first treated 
by Cashcr and Lcbowitz^^'^® and one gets J ^ j^r/j- 

(iii) In general one gets J ~ -j^ where a depends 
on the low frequency behaviour of the spectral functions 
Gpp{uj) and Gp'p'{uj')'^'^ . 

Note that the case a < 1 leads to infinite thermal 
conductivity while a > 1 gives a vanishing conductivity. 
Thus, depending on the properties of the heat baths, the 
same wire can show either superconducting or insulating 
behaviour. The usual Fourier's law would predict J ~ 
1/N, independent of reservoirs. Thus Fourier's law is not 
valid in quantum harmonic chains, even in the presence of 
disorder. This breakdown of Fourier's law in ID systems 
has been noted in a number of earlier studies on classical 
systems^^ which have looked at the effects of scattering 
both due to impurities and nonlinear ities. 



D. Temperature profiles 

The local temperature of a particle can be determined 
from its average kinetic energy , fee; = {pf/{2mi)). We 
get 

kei = - I dumiUj'^[k'^Z^l{Lo)Z[^l{-Lo)I{u;) 

+k"Z^^{u;)Z^^{-uj)l'{u;)] (25) 

This is straightforward to evaluate numerically for given 

systems and reservoirs. For the special case of heat 
transmission through a perfect one-dimensional harmonic 
chain attached to one-dimensional reservoirs through 
perfect contacts {i.e k = k' = 1) Eq. (25) simplifies (for 
large N) to 

= ^ / dkiJkifiiOk) + f'iujk)]. (26) 



where Wk = 2sin(fc/2). For T = T', we get: 



kei 



— / dfca;feCoth(-— ) 
47r ^0 2 



(27) 



which is the expected equilibrium kinetic energy density 

on an infinite chain. For weak coupling to the reservoirs, 
which can be achieved by making k and k' small, we ex- 
pect that the enc;rgy density profile for the system should 
correspond to that of a finite chain. We verify this nu- 
merically by evaluating Eq. (25) for fc = fc' = 0.1 and 
T = T' [Fig. (3)]. We compare this with the equilibrium 
kinetic energy profile of a finite chain given by: 



and the spectral function G. 



pp 



-ik 



where 



kei = i^a;,coth(^)V?(Z) 



(28) 
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FIG. 3: Kinetic energy density profile in a pure harmonic 
chain (N = 8) attached to reservoirs at equal temperatures 
T = T' = 0.2. Two different kinds of reservoirs are con- 
sidered: one-dimenionaJ reservoirs (RG) and delta^correlated 
noise reservoirs (white). The exact equilibrium density pro- 
files for an infinite chain (free) and one with fixed ends are 
also given. 



FIG. 4: Kinetic energy density profile in a pure harmonic 
chain {N = 64) attached to one dimensional reservoirs at 

temperatures T = 1.0 (left) and T' = 0.5 (right), for perfect 
and imperfect {k = k' = 0.9) contacts. The temperatures 
considered arc not very high and so the bulk temperature is 
different from the classically expected value Tav = 0.75. 



IV. DISCUSSION 



where Vs(0 = [2/(^ + l)]^^^sin{kl), ujs = 2sin(fc/2) 
where k = STr/{N + 1) with s = 1,2,...N. Note that 
unhke in the classical case where the energy density is a 
constant, in the quantum case, this is not always true. 
It is instructive to look at the equilibrium properties for 
the case where the driving is by a delta-correlated noise ( 
case (ii) discussed earlier ) . In this case the weak coupling 
limit corresponds to taking the damping constant 7 << 
1. The temperature profiles obtained from Eq. (25) for 
two different values of 7 are plotted in Fig. (3). 

We now consider temperature profiles in the nonequi- 
librium case (T ^ T'). For the Rubin-Greer ( or Lan- 
dauer case i.e. ID reservoirs, perfect contacts ), at high 
temperatures the local temperature is given by Tj = 2fce; 
and from Eq. (26) we get Ti = {T + T')/2 which is the 
classical result^^. At low temperatures and imperfect 
contacts k,k' ^ 1 we evaluate the local kinetic energy 
profile numerically using Eq. (25). As can be seen in 
Fig. 4 the temperature in the bulk still has the same con- 
stant value. At the boundaries however we see a curious 
feature noted earlier by^'*'^^: the temperature close to the 
hot end is lower than the average temperature while that 
at the colder end is higher than the average. For the case 
with delta correlated noise, at high temperatures, we re- 
cover the temperature profiles obtained ealicr for classical 
chains in^^. At low temperatures we get results similar 
to those found by Zurcher and Talkner^* and there seem 
to be some qualitative differences from the classical tem- 
perature profiles, depending on the value of 7. 



We note that the more popular approach of treating 
open quantum systems is the Caldeira-Leggett formula- 
tion. In that approach, one deals with density matrices 
and the treatment becomes complicated. In the context 
of the present problem one is not really interested in 
the full distribution but rather in physical observables 
like the steady state currents and densities and these 
are basically second moments of the distribution. The 
FKM formulation is then more appropriate and for linear 
systems one can get exact results. The other approach 
of treating nonequilibrium systems which has been used 
quite extensively in the mesoscopic context is the Keldysh 
formalism. This is a perturbative treatment where one 
writes equations of motion for a set of Green functions 
and relates them to self energies through the Dyson equa- 
tions. The current is expressed in terms of these Green 
functions. In special cases the Dyson equations can be 
solved exactly and indeed some of our results can be 
obtained'^' On the other hand our method is more 
transparent and direct. We integrate out the reservoir de- 
grees of freedom to get effective Langevin type equations 
of motion for the system. These are solved and quickly 
lead to useful results on currents and densities of both 
particle and heat which are automatically expressible in 
terms of unperturbed Green functions. The connection 
to scattering theory is also immediate and explicit. Fi- 
nally one obtains a nice physical picture of the reservoirs 
serving as effective sources of noise and dissipation. Note 
that our approach makes connections between different 
approaches such as Caldeira-Legett, Keldysh, scattering 
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theory and the transfer-Haniiltonian method. 

The FKM was earher used in studying heat transport 
in classical disordered harmonic chains and it is particu- 
larly nice that the method can be extended to the quan- 
tum mechanical regime. Earlier results on classical chains 
are then obtained as limiting cases. The more general 
quantum mechanical results can be expressed in forms 
where one can see connections with other approaches 
such as Landauer, Kcldysh, etc. 

The dependence of transport properties of a system 
on the reservoir properties is at first glance a surprsing 
fact and we briefly comment on this. From our usual 
experience in the macroscopic world, one usually thinks 
of the conductivity of a system as an intrinsic property, 
not dependent on the properties of reservoirs. Imagine 
making a measurement of the thermal conductivity of 
a wire by putting its ends in contact with heat baths 
at two diS'erent temperatures and measuring the result- 
ing current. The normal expectation is that the answer 
should not depend on the material properties of the heat 
baths. And indeed this expectation holds out true quite 
often. One physical way of understanding this is that, 
as long as the system (the wire) is a strongly interacting 
system, with good ergodicity properties, then one can ex- 
pect that, soon after contact is made with the reservoirs, 
the ends of the wire would reach a state of local thermal 
equilibrium with the reservoirs. This local equilibrium 
would be completely determined by just the temperature 
of the reservoir and this then drives the current in the 
wire. In the mesoscopic domain however there are situa- 
tions when the interactions between the carriers arc not 
strong enough to let the system reach local equilibrium. 
And then one finds that the conducting properties of a 
wire is no longer intrinsic to the wire but depends on 
details of the reservoirs. Thus any calculation of trans- 
port properties would require a detailed modelling of the 
reservoirs. An explicit demonstration, of the conditions 
under which reservoir-dependence goes away, does not 



seem to exist at present. 

As has been shown here the FKM method works as eas- 
ily for both electronic transport in disordered fermionic 
wires and thermal transport in disordered harmonic 
chains. In both cases we are able to obtain exact for- 
mal expressions for particle and thermal currents and 
these have very similar forms. Both depend on details 
of the reservoir spectral functions. The usual Landauer 
case where one writes the current in terms of transmis- 
sion factor of one-dimensional plane waves is shown to 
follow, exactly, for the choice of one-dimensional reser- 
voirs and perfect contacts. In general however one needs 
to use modified Landauer formulas and this can be quite 
crucial in interpreting experimental data. For example 
we have shown that the oscillations in conductance seen 
in the experiments by Kong et al cannot be explained 
unless the contacts and reservoirs are treated quantum 
mechanically. We also find the rather counterintuitive 
prediction that imperfect contacts can enhance the con- 
ductance of a wire. In the phonon case we make a couple 
of predictions that are interesting from the experimen- 
tal point of view: (i) the large system size behaviour of 
the heat current is a power law and the power depends 
on reservoir properties (ii) temperature profiles in per- 
fect wires show somewhat counterintuitive features close 
to contacts. It would be interesting to sec if our predic- 
tions, which are true for strictly one-dimensional chains, 
can be verified in experiments on nanowires. 
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